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Abstract. Periodic layered media can reflect strongly for all incident angles 
and polarizations in a given frequency range. Quarter-wave stacks at normal 
incidence are commonplace in the design of such omnidirectional reflectors. We 
discuss alternative design criteria to optimize these systems. 
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1. Introduction 

At optical frequencies, metallic mirrors reflect strongly for any angle of incidence 
and any polarization. However, they display dissipative losses, which constitutes a 
drawback in practical applications. Photonic crystals pQ were originally proposed by 
Yablonovitch j2j to solve this problem: they are periodically microstructured materials 
that reflect in stop (forbidden) bands within which light propagation is not possible 
in an infinite structure. Since they are made from transparent materials, photonic 
crystals can be almost free of losses at any prescribed frequency. 

In the one-dimensional photonic crystal is nothing other than a periodic 

multilayer. Although much attention has been payed to dielectric Bragg mirrors 
consisting of alternating low- and high-index layers, certain aspects of the reflection 
by periodic layered media are universal: for N periods the reflectance goes to unity 
as N~ 2 at the band edges, while tends to unity exponentially at the band gaps 00]. 

This means that, in practice, not very many periods are needed to have a stop 
band. One is then led to consider stacks of N periods (which are often called finite 
periodic structures) and apply to them conditions that are valid only when the system 
is strictly infinite. Recently we have put forward this problem and provided an 

alternative framework for dealing with these finite periodic structures: the trace of the 
basic period allows us to classify them into three classes with quite different properties. 

Concerning the performance of these structures, it is indisputable that Bragg 
quarter-wave stacks (designed for normal incidence) are the most thoroughly studied 
in connection with omnidirectional reflection [7|l%|[§i ni)ll 1 1| H 2j . In spite of this, the 
current interest in extreme ultraviolet and soft x-ray |14| optics is driving a 
great deal of work on new methods for optimizing the design of multilayer mirrors. 
In addition to the simple but cumbersome optimization by eye, only recently more 
sophisticated techniques have been started to be used |15| : relevant examples include 
the downhill simplex algorithm the systematic search in the parameters space |17| . 
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the simulating annealing |18| . the needle variation technique |19| . the Levenberg- 
Marquardt algorithm 120], or genetic algorithms 21 . The aim of this paper is to 
provide an alternative optimization criterion that has the virtue of using a simple 
analytical figure of merit with a very clear physical meaning. 



2. Notations and general relations 

We start by examining the properties of the basic period of our structure, which 
is assumed to be lossless. The field amplitudes at each side (called ambient and 
substrate) of the unit cell are related by the 2x2 complex transfer matrix M os that 
can be expressed as |2~2"ll2"3"] 

M _ ( VT as R* as /T* as \ ( a /3 \ 

Mas - \ R as /T as l/T* as ) ~ \ (3* a* J ' [L) 

where R as and T as are, respectively, the overall reflection and transmission coefficients 
for a wave incident from the ambient. Note that for identical ambient and substrate 
media we have dct M as = +1, which is equivalent to |i? as | 2 + \T as \ 2 = 1- 

We take as known the theory of reflection from multilayers 0J and its main result 
for our purposes, namely that strong reflection will occur when [Tr (M as )] 2 > 4 (these 
conditions, one for the s polarization and one for the p polarization, locate the band 
stops for each basic polarization). 

When we consider a finite periodic system that consists of N basic periods, it is 
possible to show that in the stop bands the reflectance takes the general form |4*ll51l24| 

roW = ^ (o) 

|/3| 2 + [sinh(x)/sinh(iVx)] 2 ' 1 ; 

where cosh(x) = Re(a). We are considering only positive values of Re(a) since negative 
values can be treated much in the same way. 

In practice, it is usual that the basic period of the structure consists of two thin 
homogeneous dielectric slabs with low, ni, and high, nh, indices of refraction and 
corresponding thicknesses de and dh. These Bragg structures are also appropriately 
denoted as [LH]^, where N is the total number of periods. In such a case, the condition 
[Tr (M as )] 2 > 4, required to have a stop band, can be written as 

| Re (a)\ = cos de cos dh — sin 5t sin Sh | > 1, (3) 

where Si = (27r/A)Ai is the phase shift of a wave of wavelength in vacuum A in 
traversing the layer ith and Aj is the corresponding optical path, of value 



Aj = nidi cos 9i =d i \Jn 2 i - sin 6 , (4) 

#o being the angle of incidence. For simplicity, we have assumed that the system is 
imbedded in air. 

The function Aih is 

A ih = \±^, (5) 

where r^h is the Fresnel reflection coefficient for the interface £-h. This function A«j 
is frequency independent but takes different forms for s and p polarizations. However, 
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one can check that, irrespective of the angle of incidence, the following relation for 
both basic polarizations holds: 

rrW-) 2<1 - ^ 

Att(s) \nhj 

Due to the restriction @, whenever Eq. J3J) is fulfilled for p polarization, it is always 
true also for s polarization. In consequence, the p-polarization stop bands are more 
stringent than the corresponding s-polarization ones. 

Because of historical reasons |35], typical use of dielectric mirrors has been 
evaluated at normal incidence, with layers at a quarter-wavelength thick (at the design 
frequency): 

dl = T = w 4 = T = W (7) 

The optical paths are equal and maximum reflection occurs at the frequency 

IT C TT C . . 

2 nidi 2 nhdh 
which is the center of the stop band. 



3. Optimization strategy 

In order to explain the optimization criterion we wish to propose, in Fig. 1 we have 
plotted the reflectance 9v- N > of a Bragg [LH] W structure as a function of the angle of 
incidence for several values of the number of periods N and for p polarization. The 
layer thicknesses are chosen as in Eq. Q. 

It is clear from this figure the well-known fact that the reflectance tends to unity as 
N increases. It is then hardly surprising that the usual designs found in the literature 
use the quarter-wave thicknesses at normal incidence without raising any doubt about 
such assumption. In practice not very many periods (say N ~ 10) are used in the 
visible, and one may be tempted to ask whether other thicknesses could improve the 
performance of the structure. Such a problem could be attacked by a straightforward 
computation of the reflectance as a function of layer thicknesses. With the a priori 
information that the optimum condition is close to a quarter wave stack, it will not take 
too much computational effort to find a reasonable solution [26. 27, 28, 29 . However, 
our goal here is to provide a more systematic way of dealing with this question. 

A reasonable option for optimizing the system is that, when 9$ varies from to 
7r/2, the area under the curve 1 — ISKW (which is the transmittance of the system) 
would be as small as possible. Therefore, once the materials and the wavelength are 
fixed, we propose 

«r/2 

A(d t , d h )= [1 - JftW (dt, d h , 9 )\ d9 , (9) 
Jo 

as figure of merit for the periodic structure. Alternatively, given the characteristic dip 
appearing in the reflectance for p polarization, one could also impose that this dip 
would be as smaller as possible. We have numerically checked that both criteria give 
essentially the same results. 

For s polarization there is no dip. The different behaviour at oblique incidences 
of s and p reflectances has been analyzed previously Since for a given stack the 
reflectance is greater for s polarization than for p polarization for every incidence 
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Figure 1. Reflectance for Bragg [LH]^ structures with layers at a quarter- 
wavelength thick at A = 10 (im and p polarization as a function of the angle of 
incidence. The refractive indices are ri£ = 1.75 and n h = 3.35. Each curve is 
labeled with the corresponding number of periods JV. 
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Figure 2. Reflectance for the same Bragg [LH]^ structures (with TV = 4 and 
N = 7) as in Fig. 1. The continuous lines correspond to quarter- wave thicknesses 
and the broken lines to the optimum thicknesses in Table 1. 
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angle (they are equal only at normal and grazing incidences), we argue that once the 
area A(di,dh) is optimized for p polarization, it is also improved for s polarization, 
which seems quite plausible. 

We have employed an easy-to-use quasi-Newton algorithm in order to find the 
minimum of the function (J5J), subject only to the physical conditions < dg/Xe < 1 
and < dh/Xh < 1, because the periodic character of the solutions. The numerical 
results of this optimization are shown in Table 1. Note that the use of adimensional 
thicknesses de/Xe and dh/Xh simplifies the presentation of the results, although by no 
means the results are universal: they apply only to this particular system because 
we are not considering dispersion in refractive indices (in all the paper we take the 
wavelength in vacuum A = 10/im). In any case, the proposed optimization strategy is 
independent of the particular example under study. 

For a better understanding of the performance of our optimized structure, in Fig. 2 
we have plotted the reflectance %V- N > (with N = 4 and N = 7) and p polarization for 
both the quarter- wave thicknesses as in Eq. J7J) and the optimum thicknesses in Table 
1. The improvement is remarkable: in terms of areas, we have 20 % and 100 %, 
respectively. 
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Figure 3. Contour plot of A{de, d^) as a function of the adimensional thicknesses 
di/Xi and d^/X^ for the same structure as in Fig. 2 with N = 7 periods. The 
labels in the contour lines show the value of the area (except for a factor 10 — 2 ). 
The optimum working point corresponds to the center of the grid. 

To test the tolerance of the optimum thicknesses against small deviations, in 
Fig. 3 we have a contour plot of A(de, dh) as a function of the adimensional thicknesses 
de/Xe and dh/Xh for the same structure as in Fig. 2 with N = 7 periods. The elliptical 
contours delimit the range of thicknesses giving a definite value of the area. Moreover, 
if we take the projection of the major axis of the ellipse on the coordinate axes as 
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a qualitative measure of the maximum tolerance for a given area, we conclude from 
Fig. 3 that one must be more careful in controlling the thickness dg than dh- 
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Figure 4. Contour plot of the area as a function of the refractive indices rig and 
for the same structure as in Fig. 3 but with fixed optimum thicknesses dg and 
dh given in Table 1. The labels in the contour lines show the value of the area 
(except for a factor 10 — 2 ). 

In the same spirit, one could ask about how critical is the behaviour of the area 
under deviations in the values of refractive indices. In Fig. 4 we have plotted a contour 
plot of the area as a function of ng and for the optimum thicknesses dg and dh given 
in Table 1 for N — 7. We find again the same kind of elliptical contours as in the 
previous figure, but now the variation is much smoother, indicating that, roughly 
speaking, the role of refractive indices is not so crucial as the role of thicknesses. 

Obviously, as N grows the improvement in the area is larger. In fact, tabic 
1 suggests a considerable improvement that, to some extent, may be illusory: it 
only concerns a decrease in the integrated transmission, while for the reflection, 
the improvement is not so impressive. However, the important point is that the 
improvement is just in the dip of the reflectance. We conclude finally that the method 
is especially appropriate for moderate values of N (N ~ 10), which constitute a typical 
experimental situation. 

The results presented so far hold only for the given ratio of refractive indices. To 
show that the method can be employed for arbitrary parameters of the multilayers, in 
Fig. 5 we have plotted the optimum adimensional thickness (dg/\g) op t as a function of 
the refractive indices for N — 7, while in Fig. 6 we have represented the optimum 
values of (dh/Xh) op t- We have assigned a zero thickness whenever the condition 
[Tr (M as )] 2 > 4, required to have a stop band, is not fulfilled for some values of 
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Figure 5. Optimum thickness (df/A^) op t as a function of the refractive indices 
ri£ and n/j for a [LH]^ structure with N = 7, A = 10 (im and p polarization. 




Figure 6. The same as in Fig. 5 but for the optimum thickness (rfh/A^) op t. 

the incidence angle 9q. The abrupt step is the same in both figures and gives the 
boundary of omnidirectional reflection for the stack. 

As a general feature, we note that (di/\g) op t < (dh/Xh)opt for every pair of 
allowed refractive indices. In Fig. 5, the optimum thickness varies pronouncedly with 
n£ until ne ~ 1.5 and then it is almost constant, while decreases slowly with rih in 
the region of interest. In Fig. 6 we observe that the optimum thickness decreases for 
a fixed value of ni (resp. n^) when (resp. ng) increases, and has a flat variation. 
The same qualitative behaviour is also observed for other values of N. 

In summary, we have proposed a simple figure of merit that allows for an 
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Table 1. Adimcnsional thicknesses d(/\( and d h /\ h that optimize the area for 
each value of the number of periods N. _4 op t stands for the optimum value of the 
area, while .4(1/4) represents the value obtained for quarter-wave design. The 
last column shows the percentage of improvement in these areas computed as 
100 X [.4(1/4) — .4opt]/-4opt. The data of the structure are as in Fig. 1. 



N 


de/X e 


dh/^h 


.4opt 




.4(1/4) 




Improv. 


4 


0.2567160 


0.2583615 


0.6277996 xlO" 


-i 


0.7539873 


xlO- 


-i 


20 % 


5 


0.2578014 


0.2596847 


0.2107099 xlO" 


-i 


0.2961738 


xlO- 


-l 


41 % 


6 


0.2582820 


0.2604825 


0.7195052 xlO" 


-2 


0.1199775 


xlO- 


-l 


67 % 


7 


0.2582544 


0.2610819 


0.2508925 xlO" 


-2 


0.5007563 


xlO- 


-2 


100 % 


8 


0.2578614 


0.2616020 


0.8916341 xlO" 


-3 


0.2144915 


xlO- 


-2 


141 % 


9 


0.2572066 


0.2620937 


0.3219379 xlO" 


-3 


0.9386422 


xlO- 


-3 


192 % 


10 


0.2563618 


0.2625808 


0.1177155 xlO" 


-3 


0.4179370 


xlO- 


-3 


255 % 


11 


0.2553779 


0.2630740 


0.4347006 xlO- 


-4 


0.1886905 


xlO- 


-3 


334 % 


12 


0.2542927 


0.2635773 


0.1617585 xlO" 


-4 


0.8614448 


xlO- 


-4 


433 % 


13 


0.2531341 


0.2640914 


0.6055056 xlO" 


-5 


0.3968423 


xlO- 


-4 


555 % 


14 


0.2519235 


0.2646149 


0.2277057 xlO" 


-5 


0.1841661 


xlO- 


-4 


709 % 



improvement in the performance of omnidirectional reflectors, especially when the 
number of periods is not too large. The method can be easily extended to other 
materials and wavelengths and predicts optimum thicknesses that can depart from 
the usual quarter-wave design. 
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